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1. Introduction 



Neutron scattering provides a direct measure of 
spatial and temporal correlations of atomic posi- 
tions and magnetic states of the target system. 
Thus, to a considerable extent the theory of neu- 
tron scattering is aimed at calculating physical 
quantities that are of interest beyond the scattering 
technique itself. This outline is a short and highly 
selective tour through the basic mathematics and 
concepts of the canonical neutron scattering for- 
malism, which may provide useful preparation for 
the more focused contributions that follow in this 
collection. Our approach is didactic but also takes 
much for granted. We start abruptly, for example, 
with the general formula for the differential cross 
section and along the way leave out a great deal 
that could be said; important topics are ignored, 
and some of the discussions may seem to end pre- 
maturely. If terse, however, the outline also strives 
to be logical and cohesive and to set out a di- 
gestable overview of what is, after all, a far-ranging 
discipline. 



Fortunately for the reader, excellent and well- 
known introductory articles and texts are available, 
as well as advanced treatments of both general and 
specialized interest. References [1-18] of this arti- 
cle serve as a list of some of these. Scherm's [13] 
succinct and lively tutorial on the mathematical 
formalism of neutron scattering was the inspiration 
for the present treatment, which borrows liberally 
from it. The opening chapters of Bee's book [4] 
also give a thorough and readable survey of the 
formalism. The introductory review by Price and 
Skold [10] has broad scope and provides a bal- 
anced development of both theoretical and practi- 
cal material, while the authoritative collection of 
articles edited by Skold and Price [17] reaches into 
many areas of application. Squires' clear and effi- 
cient textbook [18] is a standard and helpfully in- 
corporates mathematical preparation that other 
sources often leave out. The treatise by Marshall 
and Lovesey [8] is a classic presentation of the for- 
malism. 
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The author is indebted to colleagues who pe- 
rused drafts of this article at various stages of 
preparation. Even in its brevity it tries to convey 
their insightful comments. 

2. Basic Formulas 

2.1 The Born Approximation 

Our starting point for this outline of neutron 
scattering formalism is the definition of the double- 
differential scattering cross-section per particle, as 
obtained from the Born approximation: 

/ dV \ 1 / m y A: y , . 



" . -0)] , 



(1) 



where fl denotes the scattering solid angle, (o, the 
angular frequency corresponding to energy transfer 
h(o, m, the neutron mass and A^, the number of 
nuclei in the scattering system. Note that in some 
treatments, the differential cross-section is defined 
with respect to an energy differential, which intro- 
duces a dimensional factor of ^"' on the right-hand 
side of Eq. (1). The incident neutron beam is in the 
plane-wave state |A;o, s^, where ka is the incident 
wavevector and sq is the incident spin. The scat- 
tered beam is in the state \k,s), where k and s refer 
to the outgoing beam. The fundamental formula in 
Eq. (1) effectively is the product of a kinematical 
factor, klko, with the transition rate from Fermi's 
golden rule. The matrix element of the neutron- 
scatterer interaction V is taken in the product rep- 
resentation 



where 



\k,s,n) = \k,s)\n) 



H\n)=E„\n) 



and H is the Hamiltonian of the scatterer. The 
summation over the initial states of the scatterer is 
weighted by a probability density, p (no). The mo- 
mentum transferred to the scatterer by the neutron 
beam is ^ , where 



Q=ku-k 



(2) 



is called the scattering wavevector, or, more often 
than not, simply the scattering "vector." We 



designate the energies of an incident and scattered 
neutron as Eg and E , respectively. Then, consistent 
with the right hand side of Eq. (1), the angular 
frequency o) is determined by the neutron energy, 
Eo—E, transferred to the scatterer: 



Aco=^(kl-k^) 



(3) 



Thus positive co corresponds to excitation of the 
target, which is convenient in calculations which 
focus on the properties of the scatterer rather than 
on the analysis and configuration of the measuring 
instrument. Indeed, with respect to the probing 
neutron, w > means energy loss or "down" scat- 
tering, while ft) < represents energy gain or "up" 
scattering. 



2.2 Nuclear Scattering in Homogeneous 
Systems 

For nuclear scattering the interaction potential is 
well-approximated by the Fermi pseudopotential 



Itt^ 



V(r)=^2bj5(r-Rj) , 



m 



(4) 



where bj is the isotope- and spin-dependent scat- 
tering length of the neutron-nuclear interaction, 
and the summation is over all of the nuclei in the 
system. The scattering lengths in Eq. (4) are de- 
fined for fixed nuclei in the laboratory frame, and 
are generally referred to as "bound " (or "bound- 
atom") scattering lengths. For unbound nuclei, the 
bound scattering lengths are replaced by "free" 
scattering lengths. 



fl/ = 



Ai + 1 



bi 



(5) 



where the Ai are the nuclear atomic weights. The 
prefactor in Eq. (5) stems from the reduced mass 
of the interacting neutron-nucleus pair. The bound 
scattering lengths, in fact, are always the correct 
choice, provided the scattering is calculated cor- 
rectly, as illustrated below by the case of the ideal 
gas. Thus tabulations of scattering lengths (e.g., the 
comprehensive table by Sears [14]) give the bound 
values. For most nuclei, of course, the differences 
between the bound and free scattering lengths are 
relatively small, but hydrogen and deuterium are 
notable exceptions. 
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For a neutron beam of arbitrary spin polariza- 
tion, the scattered intensity is determined by 



_d^ 
d/3 



£=2:p(^o)(d7^)^^^^ (6) 



where /j(i)-/7(-i) is the polarization state of the 
beam. For an unpolarized beam, /j(±i) = i, while 
for nonmagnetic interactions the cross-section is 
independent of So . When there is no knowledge 
connecting the internal (isotopic and spin) states 
of the scattering nuclei to their positions, the 
observed cross-section is an average over the 
nuclear degrees of freedom: 



Equation (8) can be rewritten even more compactly 
as 



bibj = 2 ^fjbl 



(13) 



where 



Uj = 



1 fora = "coh", 
8,j fora = "inc". 



(14) 



The combination of Eqs. (1), (4), (7), (13) and (14) 
then gives 



d^o- _/ dV \ 
d/2da> \d/2d<u/„uc 



(V) 



d'o- k 



d/2d<u 4'n-ko 



2(r''Sa(Q,(o), (15) 



It is generally an excellent approximation to 
assume that there are no correlations between the 
internal states of different nuclei. As a conse- 
quence, 

bibj = bibj, iT^j, 

where the overbar denotes the average over 
nuclear internal degrees of freedom. Then 



bibj =b?„c^ij+b^h , 



(8) 



where 



where 



•1 

Sa(Q ,to) =-r^ 2 P («o)2 ?U 



N 



n o,n I, J 



(«o|e-'e*^|«)<«|e'e«'|«o>8(^^-^ -w) . 



Explicitly, these are: 



(16) 



bcoh = b, and 
bL = P-b^ , 



(9) 



define the interactions for nuclear coherent and 
incoherent scattering. The effective scattering lengths, 
biac and bcoh , are often expressed in terms of bound 
cross-sections as. 



where 



cTa = 4'n-ba , 



•coh" , 



(10) 



(11) 



5coh(fi, to) = Tp 2/'(no)2 



N 



n o,n I, J 



(«fl|e-'e*'|«)(«|e'e«'|«o)8(^^^-^ -<o) 



and 



5i„c(fi,(o)=T7 2/'(no)2 



N 



(17) 



Therefore, 



Criac + CrcoW = 4TTb^ . 



(12) 



Kn|e'e-^K)P8(^^^-.,) . (18) 
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Notice that the site summation in 5coh (S,w) is 
unrestricted and thus includes the summation in 
5inc(S,w). The total cross-section is obtained from 



O'cotal 



=ir 



d^ 



d/2da) 



d/2d« 



(19) 



In the upper limit, Eo = ^^ld/2m is the incident 
energy, the largest possible neutron energy loss. 
Often, however, w-integral is well-approximated by 
taking the upper limit to « . The quantity 



■i — 7^ O'total , 



(20) 



for A'^ scattering nuclei in sample volume V, is 
called the macroscopic cross-section and is just the 
inverse of the scattering mean free path. 

The total cross-section, Eq. (19), depends on the 
incident n eutron energy; note, for example, the fac- 
tor ko = \'2mEo/^ in Eq. (1). The £o-dependence 
of total cross-section can be obtained analytically 
for the case of neutron scattering by an ideal gas. 
The result shows that o-totai « jBo""^ for Eo^*0, as 
expected from the kinematics, while for Eo > kaT/A , 
the asymptotic scattering is 



fcotal 



- rTTT I (''''■"= "^ '''"'h) > 



which is just the free cross-section. This renormal- 
ization of the bound cross-sections at high incident 
energy is brought about by the recoil kinetic energy 
of the impulsive neutron-atom collisions and is 
observed in condensed matter. 

The relation in Eq. (15) is the basic law of non- 
magnetic neutron scattering from homogeneous 
systems. Generalizations to heterogeneous nuclear 
scattering and magnetic scattering, described 
below, depart only slightly from this archetype. 
Thus (nuclear) neutron scattering is a sum of 
coherent and incoherent contributions, each of 
which conveys different kinds of information about 
the scattering system. Only Scoh(Q,(o) expressly 
depends on the relative positions and motions of 
the nuclei and thus provides explicit knowledge of 
structure and collective dynamics. In contrast, 
5inc(fi,w) discloses no spatial structure and only 
reveals the motions of individual particles. For 
example, while coherent scattering can measure a 



phonon dispersion curve, w vs g , incoherent scat- 
tering can only measure a phonon density of states, 
N(o)). Needless to say, a density of states is often 
just the information desired. 

The weights of coherent and incoherent scatter- 
ing in the measurement are determined— for a 
homogeneous system— by o-coh and o-jnc. which vary 
randomly from one nuclear species to another. 
Most nuclei are stronger coherent scatterers than 
incoherent scatters, and a few even scatter only 
coherently. For example, if we define 



A = 



ftotal 



(21) 



then/coh=l for the abundant isotopes of Al, C, and 
O, while /coh^O.97 for Mg and Fe. Some nuclei 
scatter mostly incoherently. For example, 
/inc(V) = 0.995 and /i„c(Co) = 0.84. Hydrogen, 
/inc(H) = 0.98, is a singular example of an incoherent 
scatterer. Its total cross-section, o-totai(H) = 82 bams 
(1 barn = 10"^m^), is by far the largest among the 
elements and, of course, it is abundant in many 
materials. Hydrogen in hydrated compounds and in 
hydrocarbons causes significant incoherent scatter- 
ing, which often masks the desired coherent scatter- 
ing in structure determinafions. Indeed, in small 
angle neutron scattering, water scattering, which is 
virtually independent of Q in the range of interest, 
is used to calibrate area detectors for spatial 
nonuniformity. Deuterium scatters incoherently, as 
well, but actually is a stronger coherent scatterer, 
with /coh(iD) = 0.74. Hydrogen and deuterium also 
are readily dissolved in many metals as a highly 
mobile interstitial impurity. Incoherent neutron 
scattering thus probes hydrogen vibrations and 
transport in these metals: in Siac(Q,(o), the 
posidon operator, Ri (t), then refers to only to H if 
the host metal is a coherent scatterer. 

In heterogeneous nuclear systems, Eq. (15) is 
replaced by 

^=|i:5r(e,.), (22) 



where 



5r(e,co)=ii:/;(«o)2H?y 

n o,n ij 

(«o|e-g*/|n)(n |e'g'*'|no)8( ^'' ~^"" -co) , 



(23) 
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and where the ^-function defined in Eq. ( 14) has 
been extended to the H-function, 



E^-= ■ 



fc«.h,/6cohj fora = "coh", 






for a = "inc" . 



(24) 



Thus, in heterogeneous systems, including crystals 
with bases, the coherent scattering is also sensitive 
to the signs of the constituent fecoh-values, which 
also vary randomly with species, although bmb > in 
most cases. For example fecoh(H) = — 0.37, while 
6coh(D) = 0.65. This sign difference is exploited in a 
wide variety of neutron scattering experiments, in 
which deuteration is used to label molecular sites 
or adjust overall coherent scattering contrasts to 
enhance the contributions from various compo- 
nents of the structure. 



2.3 Van Hove Formalism 

The Van Hove formalism recasts the Bom 
approximation formula for Sa(Q,ci)) into the 
Heisenberg representation and introduces quan- 
tum generalizations of classical correlation func- 
tions. Begin with the identity. 



y- [ e'"' {m\Aop(t)\n)dt , (25) 



where 



H H 



■^op\j ) — e A A ope ft 



(26) 



Then Eq. (16) becomes 



5„(fi,6>)=^J ^e-'"7„(e,0d/ , (27) 



with the introduction of 






(28) 



Two additional steps were also needed to reach 
Eq. (28). First 



e'^e''C*e-'^ = e-''fi'<') 



and 



5]p(«o)<«oMop|no> = <^op> . 

"0 

Finally, with the use of the obvious identity. 



jQRs = j e'e'-'5(r-J?;)dV , 



Eq. (28) takes the form 



I.{Q,t) = \d^-'>G.(r,t)d'r . 



(29) 



where 



{d(r' -Ri(0))d(r' +r-Rj(t)))d'r' . (30) 



Then Sa(Q,<o) and Ga(r,t) are connected by the 
Fourier transform 



5.(6,..)=^ J £ 



e'(fi'-"')G„(r,r)dVdf , 



(31) 



and its inverse, 
Ga(r,0= J J_ 



e-'«-'-)5„(fi,6,)-^d6, 



(32) 



The asymmetric placement of the lir in these 
Fourier integrals is unfortunate, perhaps, but con- 
ventional. It could be avoided, say, by multiplying 
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the right-hand-side of Eq. (28) by Itt, but this is 
never done. The relationships among the quantities 

Sa(Q,(o)-^ "scattering function", 

Ia{Q,t) -* "intermediate scattering function", 

Ga{r,t) -^ "correlation function", 

comprise the Van Hove representation of the Born 
approximation. Now the nuclear position opera- 
tors, Ri(t), are not constants of motion, so that 

[R>(t),Rj(t')]^0 , 

and the substitutions implied by the the 8-functions 
can not be made in Eq. (30), except in a classical 
approximation. However, the correlation function 
in the coherent case can be exactly represented by 
the form 

Gcoh(r,0=^J(p(r',0)p(r'+r,0)dV' , (33) 

where 

p(r,t) = 2S(r-R<(t)) (34) 

i 

is the nuclear particle density operator. Also, with 

p(G.O = Je'^Xr,OdV, (35) 

the coherent intermediate scattering function is 
also represented by 



/coh(G,0 = ^(p(-G,0)p(e.O)- 



A^ 



(36) 



In materials applications, the meaning of ( • • • ) is 
generalized to include configurational averaging, as 
well as the implied thermal average. For a heteroge- 
neous scattering population, bcobp(r,t) becomes 
the scattering-length density field, 



V(r ,t)= ^ bcoKiSir -Ri(0) , 



(37) 



and Scoh(Q,<o) is a sum of contributions depending 
on self- and cross-correlations among the various 
nuclear populations. 

The Van Hove formulation was a signal advance 
in non-relativistic scattering theory and a major 
stimulus to condensed matter applications beyond 
crystallography. Although rooted in Fermi's golden 
rule, it manages, with the introduction of the corre- 



lation function, to "hide" its overt quantum 
mechanical origins and bring forth a classical 
looking representation that is a powerful guide to 
intuition, while also providing the path to relation- 
ships with statistical mechanics and linear response 
theory. 

In thermal equilibrium, the statistical weight 
p (no) of the initial states is the Boltzmann factor 



-/5£nO 



where 



/?(«o)=— 2" 



-pEn 



The energy conservation explicit in Eq. (1) means 
that 

p{no)-^p{n)c^''°'. 

within the sums over states. This leads eventually to 
the principal of detailed balance : 

S.{-Q,-(o) = &-^"-'Sa{Q,o>). (38) 

In applications to inelastic scattering, Eq. (38) pre- 
dicts that "down" scattering (w>0) is always 
stronger than "up" scattering (<t)<0). For centro- 
symmetric systems, the scattering function is an 
even function of Q , 



Sa(-Q,co) = Sa(Q,<o). 



(39) 



Actually, this is also true for systems without a cen- 
ter of symmetry if the inverted equililbrium struc- 
ture differs from the original by a constant 
displacement. 

The scattering function, Sa(Q,<o), is defined to 
be real, but the derived correlation function, 
Ga(r,t), is complex: 

G(r,t) = G'(r,t) + iG"(r,t). 

This can be seen from the integral in Eq. (32), since 
Sa(Q,(o) is not an even function of w— the con- 
straint of detailed balance. Thus a general relation- 
ship exists between the real and imaginary parts of 
the correlation function, which ultimately leads to: 



2i 



I gCe— ')G;'(r,OdVdf . 



(40) 
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This leads in turn to a concise connection between 
the scattering theory and linear response theory: 

S^.(Q,<o)= ~^_^,,^ x"(Q,<o), (41) 

1 - e '^ 

where xiQy^) is the generalized susceptibility 
function for the density operator in Eq. (34). 

2.4 Magnetic Scattering of Unpolarized 
Neutrons 

Since the neutron has a magnetic moment, inci- 
dent neutrons interact not only with the nuclei of 
the atoms in the scatterer but also with the mag- 
netic moments of their unpaired electrons. For un- 
polarized beams, the magnetic interaction with the 
nuclear spins is already incorporated into the dis- 
tinction between bmii and biac in Eq. (9). The poten- 
tial for the neutron-electron interaction is 



F™g(/-)=-M-fi(r), 



(42) 



where fi = jfiNtf is the neutron magnetic moment 
operator, in terms of the neutron gyromagnetic 
ration, 7 = 1.9132, the nuclear magneton, /an, and 
the Pauli operator, a. The internal magnetic field 
operator, B(r), is 

B(r) = 2^.2[^x{^)+^], (43) 

where fis is the Bohr magneton, s, and p, are the 
electron spin and momentum operators, respec- 
tively, and the summation is over all unpaired elec- 
trons. The (") symbol denotes a unit vector. For the 
Born approximation the required matrix element 
offi(r) is 

x-,r- •- i(PiXQ)'] 

{k \B (r)\ko) = 2mb2 [e X (s, X Q) + \q \ ^'^' 

(44) 



= AioM(e), 



(45) 



where /ao is the vacuum magnetic permeablity and 
M is the magnetization operator. For an unpolar- 
ized beam, the differential cross-section can be 
expressed as in Eq. (15) with the identifier a taking 
on the additional value "mag", and with the defini- 
tion 

o-mag = 4Tr(y/-o)^ (46) 

where /"o is the classical electron radius. Note that 
this usage of the symbol, o-mag, is not universal. The 



generalization of the Van Hove representation to 
the magnetic case gives. 



5mag(e,w) = 2 



S^v — QfiQv 1 



f I^^'{Q,t)dt, (47) 

J — 00 

where /a and v are coordinate labels. The interme- 
diate scattering function I^v^(Q,t) extends the 
definition in Eq. (36), 

/™^ (Q,t)~{M^(-Q,0)M.(Q,t)) , (48) 

for //m magnetic ions, while the corresponding mag- 
netization correlation function, from Eq. (33), is 

G;r.^«(r,O=;^j<M,(r',0)M„(r'+r,O)dV'. (49) 

Thus, except for the appearance of the quantity 
known as the Halpern tensor. 



O^u i^t^l^v , 



(50) 



in Eq. (47), the formal extension of nuclear scatter- 
ing to magnetic scattering of unpolarized beams 
simply entails replacing the scalar nuclear density, 
p, by the components of the magnetization, M,,,. 
The Halpern factor, however, is responsible for the 
special, anisotropic character of magnetic neutron 
scattering. For example, for a single-domain ferro- 
magnetic sample. 



8^v-Q^Qv-^l-(Q-My 



(51) 



where M is the direction of the domain magnetiza- 
tion. Thus, in the Born approximation, there is no 
magnetic scattering for g || M . As indicated by the 
scattering law, Eq. (15), now with 



a= ' 



'coh", 
'inc", 
'mag". 



(52) 



nuclear and magnetic scattering coexist, in general, 
in magnetic materials and can have comparable 
weight. The anisotropy of magnetic scattering offers 
one means of separating these contributions by 
choosing directions of Q along which the magnetic 
component is suppressed. Squires [18] introduces 
basic elements of magnetic scattering, while Balcar 
and Lovesey [3] provide a thorough exposition. 
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3. Generic Applications 
3.1 Elastic Scattering 

The scattering function can be analyzed into 
distinct w-dependent contributions according to: 

5„(G,co) = /„(G,»)5(o)+5r(G,co). (53) 

The first part is called the elastic scattering, which, 
by convention, means sQattering with no change in 
the neutron energy, so that k =ko. The isolation of 
the elastic contribution is the result of writing 

in the integrand of Eq. (27) and assuming that the 
limit, 

/„(G,») = Hm[/„(e,0], 

exists. The first term in Eq. (54) is independent of t 
and thus produces the scattering proportional to 
S(a)). The «o«-elastic part SS'(Q,ot)) stems from 
the Fourier transform of the r -dependent remain- 
der and comprises both /« elastic and quasie\astic 
contributions. Furthermore, for |r-/'|-» <», 



(A(t)B(t')} = {A(t)}{B(t')} 



(55) 



where the A "> is equilibrium average of the opera- 
tor/!. Therefore, in Eq. (53) 



PMr)= J{5(r')){5(r'+r))dV' (58) 

can be interpreted as the conditional probability 
density for nuclear position r at time t—*oo, given 
that r = at r = 0. These expressions for the elastic 
scattering are easily extended to heterogeneous 
systems. 



3.2 Energy-Integrated Scattering: The Static 
Approximation 

Coherent elastic scattering, which reveals equi- 
librium structure, must be measured by energy- 
resolved techniques. Often, however, structure is 
studied using methods that do not discriminate the 
scattered neutron energy; these experiments 
observe the energy-integrated cross-section 



do- f 
d/2 J. 



dV 
dfidco 



dct> . 



(59) 



In the static approximation , this is represented as 

where 

Sa(Q)=\ S.(Q,(o)d(o. (61) 

"' — CO 



/coh(e,«)=^|p"'(e)r. 



(56) 



which is the basis of structure determinations in 
solids. Various aspects of this extensive field are 
elucidated in [1], [7], [8], [11], and [18]. For liquids, 
on the other hand, p'\r)=NIV, so that 



K(e,»)=^(p)'5(G), 



leading to the dictum that liquids scatter elastically 
only in the forward direction. The weight of inco- 
herent elastic scattering is given by 



/mc(G.») = |e'°'/'scif(r,»)dr 



(57) 



where 



Since Q is held fixed in the w -integration, this 
approximation entails neglecting the G -dependence 
of the final neutron wavenumber, 

A: = |A:o-G|. 

and setting k =ko, as in elastic scattering. Then 
from Eq. (27), 

S.(Q) = h(Q,0). (62) 

Since [Ri(t), Rj(t)] = 0, this can be analyzed as 



'5a(G) = Je-'e-^G„(r ,0)d^ 



(63) 



where 



G„(r,0) = ^ S ^r.i {d{r+R>{0)-Rj{m- (64) 

i.i 
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Therefore, for the incoherent contribution, 

GiBc(r,0) = S(r), (65) 



so that 



5i„c(e)=i 



(66) 



For the more interesting coherent component. 



where 



5cx>h(e) = S(r) + g(r), (67) 



Bir) = l{d(r+Ro(0)-Rjm) (68) 

is the static pair-correlation function. Thus energy- 
integrated coherent scattering measures an aver- 
age of instantaneous "snapshots" of the time- 
dependent structure, while elastic coherent scatter- 
ing measures the equilibrium averaged structure. In 
harmonic crystalline solids, where the atoms more- 
or-less "stay put," both methods yield the equi- 
librium positions and mean-square thermal 
displacements, but the two forms of measurement 
differ in important ways. Indeed, from Eq. (53), one 
has 



5coh(G) = 5Sh(G) + 5c"o1,(G), 
where we have defined 

5^h(e)=/ooh(fi,"), 

and 



(69) 



(70) 



55fh(e)=f 5c1fh(e,<")dft>. (71) 

J — 00 



For monatomic harmonic crystals, in particular, 



(27r) 



sMQ)=-ji^^-'"'lHQ-K), 



(72) 



where Hq is the volume of the unit cell, the K are 
wavevectors of the reciprocal lattice, and 
exp(-2W') is the Debye-Waller factor, which 
depends on the thermal displacements and is 
defined below at Eq. (107). The S-functions in 
Eq. (72) are the Bragg peaks from which the equi- 
librium structure is deduced. In this context the 
integrated nonelastic scattering in Eq. (69) is 
referred to as thermal diffuse scattering. For the 



special case of an Einstein solid (dispersionless 
vibrations), the thermal diffuse scattering varies 
monotonically with Q. In general, however, for 
dispersive vibrations (acoustic phonons), the 
thermal diffuse scattering is a structured back- 
ground, including temperature-dependent "wings" 
centered on the Bragg peaks. Thermal difftise 
scattering is discussed again in Sec. 3.5.2. 

In strongly anharmonic solids and systems where 
structure changes on the time scale of the measure- 
ments, elastic and energy-integrated scattering 
probe substantially different correlations, and only 
the latter provides the short-time information 
usually of interest in such cases. 

3.3 Sum Rules 

Taken in a larger context, the formula in Eq. (61) 
is the zeroth-order member of a set of frequency 
moments of the scattering function, Sa{Q,(o), 
which in normalized form is defined by 



{co"UQ) = 



In particular. 



1 



5»(e) 



r 



o>"Sa(Q,o>)do>. (73) 



{<o'UQ) = i. 



(74) 



While seeming trivial, in practice Eq. (74) actually 
constrains physically distinct measurements and 
thus provides a means of checking their mutual con- 
sistency. The higher-order rules are formally gener- 
ated starting with Eq. (27): 



£-JaiQ,t)\.^o = i" J Ja(.Q,co)co" dco . (75) 
Then it is easy to get 

NSaiQ) ij 

^e-'e-*.(«)^e'e«^(") , (76) 



where, following Eq. (26), the derivatives with 
respect to t are determined by 
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le'e-*'W = |[//,e'e-*^(')]. (77) 



dt 



For velocity-independent interactions, the commu- 
tator is independent of the atom-atom interaction 
potential and gives 



[//,e'e-*>(')] = [-^,e'e*>(')] = 

e.fi-/<.w(^+^e.ij.(,)), (78) 



where M is the atomic mass. The contribution from 
the second term averages to zero in Eq. (76), and, 
recalling the normalization in Eq. (73), it follows 
easily that 






(79) 



which can be interpreted as the average recoil 
kinetic energy at fixed scattering vector Q. Note 
that if the f -derivative is treated classically, the 
result is equivalent to keeping only the second term 
in Eq. (78), which gives (w)a(2) = 0. This is incon- 
sistent with the asymmetry of Sa(jQ,(o) required by 
detailed balance, Eq. (38), except in the high- 
temperature limit. 

Higher moments are increasingly system-depen- 
dent and more difficult to calculate. For a classical 
liquid, the second moments are 



(■"''"(e)-^ 



The coherent case. 



{<o'U(Q) = 






(80) 



(81) 



and its analogues come into play in the analysis of 
quasielastic scattering (Sec. 3.6) in liquids and 
lattice fluids, where (w^)coh(2) is used as a measure 
of the width of the coherent contribution to the 
quasielastic line. Then Eq. (81) predicts that the 
width has a minimum near the maximum of 
Scob(Q), the effect known as de Gennes narrowing, 
which physically expresses the relative stability of 
density fluctuations that are consonant with the 
dominant short-ranged liquid order at Qmax- 



3.4 Small Angle Neutron Scattering (SANS) 

Small angle scattering is the discipline of study- 
ing microstructure of materials by investigating 
scattering in the immediate neighborhood of the 
incident beam. For neutrons with wavelengths 
larger than the Bragg cutoff —\.&., Ao>a, where a is 
the largest unit cell dimension— Bragg scattering 
from the atomic-scale structure is suppressed, and 
all scattering is near the forward direction. The 
(2 -dependence of small angle scattering thus arises 
from structural variations over length scales usually 
much larger than the atomic scale. An informal but 
helpful rule is that microstructure on length scale L 
produces scattering in the Q -range 



e<2TrL-'. 



(82) 



In other words, the more homogeneous the 
microstructure, the smaller the Q -range in which it 
is observable. To a good approximation in SANS, 
the scattering angle is given by 2d = \oQ/2TT. 
Combining this with the criterion in Eq. (82), gives 



29^ 



L • 



(83) 



This shows that increasing the neutron wavelength 
magnifies the observable scale of microstructure; 
i.e., at a given scattering angle (position on the 
detector relative to the incident beam), larger 
homogeneous regions contribute to the scattering at 
larger incident wavelengths. Small angle scattering 
also occurs for neutron wavelengths below the 
Bragg cutoff, but then intensity at small angles is 
reduced by Bragg scattering into large angles. In 
thick specimens, especially, multiple Bragg scatter- 
ing can contaminate scattering in the small angle 
regime, a bane of the early history of the field. 

As with other methods of structure determina- 
tion, small angle scattering can be implemented as 
elastic scattering, following Eqs. (53) and (56), or as 
energy-integrated scattering, as defined in the 
preceding section. Most SANS instruments, in fact, 
integrate the scattered energy transfer, and the 
static approximation usually is appropriate. Thus, 
for heterogeneous systems Eqs. (84) and (63) are 
generalized as 



(S..,=2«"(G), (84) 



and 
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5^'(e)= f e'«-'Gr (r,0)dV 



(85) 



Since the notation is getting cumbersome, for the 
remainder of this section the superscript "het" will 
be omitted, and the specification, ? = 0, will 
become implicit. Also, recalling Eq. (9), beoh will 
be denoted by the more compact b . Now, in order 
to make contact with conventional SANS terminol- 
ogy, introduce the function r(r). 



Then, from Eqs. (33) and (37), 



(86) 



r(r)=^j{r,(r')r,(r'+r))d'r'. (87) 

For the incoherent contribution, Eq. (30) easily 
gives 

Gi„c(r)={bL)S(r), (88) 



where 



(feinc)= wS&in 



N' 



(89) 



This is just the extension of Eq. (65) to hetero- 
geneous systems. Small angle scattering usually is 
analyzed in terms of the Porod-Debye correlation 
function, as defined by 



y(r) = 



r(0)-r(oo) 



(90) 



By construction, it follows immediately that 

■y(O) = 1, and 
y(oo) = 0. 

Using Eq. (87), the limits of y(r) are 

y(0) = {V%, and 



(91) 



where 



y(00) = (77)2^, 



(•••>K = ^J(-'-)dV, 



(92) 



(93) 



means a combined thermal-configurational-volume 
average. In most theoretical SANS treatments these 
averages are treated as being redundant, and the 
choice of which one to do "first" is usually a matter 



of mathematical convenience. The Porod-Debye 
correlation function can then be written as 

i/(AT7(r')AT7(r'+r))dV' 

y(r)= ^t-t;^ , (94) 



where 



{(^V)% 



^V(r) = v(r)-{v)y, 



(95) 



is the fluctuating part of the scattering-length den- 
sity field. The combination of these results gives 

N{^)^^^^^^ = N{bL) + (2^rV{(Av)%S(Q) 

+ F((A77)\Je'fi-'y(r)dV. (96) 

The first term is a (2 -independent background, and 
the second is the forward coherent scattering, 
which usually is masked by a beam stop. The 
Fourier transform of y(r) in the last term de- 
scribes the scattering caused by the spatial varia- 
tion of the microstructure. Equation (96) is the 
basic law of small angle neutron scattering. 

A typical class of materials in SANS studies 
consists of two or more homogeneous phases sepa- 
rated by sharp interfaces and partitioned into 
multiply-connected volume compartments. For 
such systems the scattering-length density can be 
represented as 



where 



17 (/•) = 2 b„p„ (r) , 



Pn(r) = 2s(r-Ri), 



(97) 



(98) 



and b„ is its coherent scattering-length, which can 
be computed as a cellular average for polyatomic 
phases. The component densities satisfy 



\Pn)y=y<Pn 



(99) 



where </>„ is the volume fraction of the n -th phase, 
and also have the properties, 



I yPn{r) 



p„(r)p„,(r) = \ iv c (100) 
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The first expresses the spatial partitioning of the 
phases, while the second is the normalized idem- 
potency condition. It then follows that 

{i^v)\ = {^)\{bl^)-{b.o.f], (101) 



where 



bcah = 2j <i>nbn 



(102) 



Then, using the phonon representation of the 
displacement operator, one can derive (with some 
tedium) the basic identity of the harmonic approxi- 
mation. 



g-e.{«ffl-«/'.)g-2H'g2H'y(g^^) 



(105) 



For the special but important case of two phases, 
this reduces to the symmetrical formula 



where 



((At, )\ = 0(1-0) (7,1 -7,2)', (103) 



w^;(G,O=^<fi-«.(0)G-«/(O>, (106) 



where is the volume fraction of either one of the 
phases, and t,i and 172 are the scattering-length 
densities. 

3.5 Inelastic Scattering 

In Eq. (53) the non-elastic contribution, 
ST{Q,o}), may describe /nelastic scattering, 
5«a«elastic scattering. Sec. 3.6, or both, depending 
on the system being studied. Inelastic scattering 
means energy-resolved scattering not centered on 
w = 0. Coherent inelastic scattering typically con- 
sists of well-defined lines riding on a diffuse, or 
structureless, ^-dependent background. The lines 
result from resonant scattering of the neutron by 
long-lived elementary excitations of the system, 
such as phonons and magnons (in magnetic materi- 
als) and impurity atom vibrations in host-guest 
systems. Diffuse coherent inelastic scattering is 
produced by scattering from multiple, or non- 
elementary, excitations. Energy-integrated inelastic 
scattering is usually the main source of the thermal 
diffuse background in structure determinations 
using energy-insensitive detection methods. Inco- 
herent inelastic scattering measures excitation 
densities of states, to a first approximation. 

3.5.1 Scattering by Phonons 

Scattering by phonons provides the prototype for 
neutron inelastic scattering. Squires [18] gives a 
lucid treatment. To calculate the intermediate scat- 
tering function, Eq. (27), resolve the nuclear posi- 
tion operator into the equilibrium coordinate and a 
time-dependent small displacement, 



and 



W = W.i{Q,Q) 



(107) 



R,{t)=Rp + Uiit)^ 



(104) 



The ubiquitous exponential, t\p{ — 2W), is the 
Debye-Waller factor. It is a long-standing conven- 
tion to leave the g -dependence of the W function 
implicit. A concrete formula for W follows from 
Eq. (113) in the next section. Although correct 
quantum mechanically, the result in Eq. (105) actu- 
ally can be obtained more directly by treating 
Q • Uj (t) as a classical Gaussian random variable, C, 
with zero mean, and by using the concommitant 
rule: 

(eO = e"^<^>. 

The Taylor expansion of the final exponential in 
Eq. (105) leads to a systematic development of the 
inelastic scattering function, 5i"°'(g,w). The first 
term (unity) gives the elastic scattering from the 
equilibrium structure, as in Eq. (56). The next 
term— the first f-dependent one— describes one- 
phonon scattering, while the remainder produces 
multi-phonon scattering. For a Bravais lattice the 
result for coherent one-phonon scattering turns out 
to be: 



,._p,.„ . (2^)^;^e-^'^ [g-6,(g)p 
{[ns{(i) + l-\S{o>-o>,{(l))S{Q-q -K) 



+ ins{q)'\8{o} + 0}s{q))8{Q -^q ■\-K)} . (108) 
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In this formula M is the nuclear mass, /2o is the 
volume of the lattice unit cell, and Kis a reciprocal 
lattice wavevector (often denoted as Ittt, where t 
is a reciprocal lattice vector). The phonon modes 
of wavevector q and polarization s have frequencies 
(Os(q), polarization vectors Ci(g) and thermal 
occupations ns{q). The two terms in Eq. (108) 
describe, respectively, scattering with the emission 
of a phonon— the neutron is down -scattered by 
energy ^(Os{q)— and scattering with the absorption 
of a phonon— the neutron is wp -scattered by 
energy fi(>h{q)- In each case the scattering satisfies 
the kinematical constraints 



£^{kl-k^)^±fico.{q) 
and 



(109) 



application of inelastic scattering. For example, in 
hydrogen in metals, the point symmetry of occu- 
pied interstitial sites usually is revealed simply by 
assignment of the intensities produced by vibra- 
tional modes having signature degeneracies. 

3.5.2 Thermal Diffuse Scattering 

With the material now at hand, we continue the 
description of thermal diffuse scattering begun in 
Sec. 3.2. In general the energy-integrated non- 
elastic scattering that constitutes thermal diffuse 
scattering incorporates both quasielastic and in- 
elastic scattering. Here, however, we only consider 
the inelastic contribution from phonon scattering. 
Thus using Eq. (105) in Eqs. (28) and (62), and 
recalling Eq. (71), one gets 



k^-k=K±Q. 



(110) 



For incoherent one-phonon scattering, the scat- 
tering function can be obtained from Eq. (108) by 
making the formal replacement. 



(27ry 
/2o 



d{QTq-Ky^ 1. 



This leads to 

5ij,c'"'(e.ft') = 

3JV^e-^"' {[Q'^s{q)?ns{q)). 



IM 



(O 



•^(<o), (111) 



where ^"{(0) is the one-phonon density of states 
and ( • • • )a, denotes the average over the surface in 
reciprocal space on which (Os(q) = (o. 

For multi-phonon processes the corresonding 
kinematical relations are not mode-specific, and 
diffuse w-dependent scattering generally is the 
result. Often the multi-phonon background must 
be estimated in order to achieve a reliable analysis 
of one-phonon scattering. Occasionally multi- 
phonon scattering may even appear in the guise of 
"broadened" one-phonon scattering in compli- 
cated spectra. 

Inelastic scattering is also used to measure spin- 
wave spectra in magnetic materials. The kinemati- 
cal constraints for 1-magnon scattering are the 
same as in Eq. (110). Determination of local vibra- 
tions in host-guest systems is another important 



55fh(e) = 1-e-^"'-^^ 



-2W 



N 



y g-;e-(«,M -«,«) (g^wvte.o) _i\ 



(112) 



where Wi,(jQ,t) was defined in Eq. (106). Explic- 
itly, forr=0. 



'j\^^ / 4MN "^ o)s(q) 

q,s ■• VI / 

\lns {q) + l]cos[? • {Rr -/?/" )] . (113) 



The formula for W, Eq. (107), is obtained by set- 
ting i =i here. For an Einstein model the thermal 
displacements on different sites are uncorrelated, 
so that Wij{Q,Q) = Q for i-^j . In this case, there- 
fore, the term in Eq. (112) with the sum over sites 
vanishes identically, and the thermal diffuse back- 
ground is l-exp(-2f^). Notice that in the opposite 
extreme, if Wij (g ,0) were independent of the dis- 
tance between sites, the summation in Eq. (112) 
would produce Bragg peaks, as in 5mh(6)- Instead, 
Wii{Q,Qi)-^Q as \Rn-Rp\^oo, and this fall-off 
moderates the site summation. However, Wij{Q,Q) 
decays slowly with distance, which leads to a weak- 
ened expression of the reciprocal lattice structure 
in the thermal diffuse scattering. This is seen by 
expanding exp[2W'y (g,0)], in analogy to the 
phonon expansion of S'coh{Q,co). Indeed, one can 
obtain the 1 -phonon thermal diffuse scattering 
directly by integrating Eq. (108) with respect to oi. 
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This gives 



SLl'\Q) = 






where 



- 2 [2n.(g) + l] 
f«(9)] 7^7 — 



q = Q-K 



(114) 



for each K in the sum. Thus the acoustic phonon 
branches give rise to algebraic singularities coinci- 
dent with the Bragg peaks, since (Os{q) = Csq for 
small \q\. Specifically, for g^O near a particular /if 
and at temperatures such that fi(D{q)<k^T, then 



SU\Q) '=^\Q-K\ 



(115) 



since ns{q)~k^T I fi(x)s{q). On the other hand, if 
^0}s(q)>kBT, then 



Siol^\Q) a \Q-K\' 



(116) 



Although these results are derived in the harmonic 
approximation, they correctly imply that soft 
modes in plastic crystals also produce significant 
thermal diffuse scattering near Q=K+qc, where 
IDs (qc) = at the softening temperature. 

Historically, before the ascendancy of inelastic 
neutron scattering, x-ray thermal diffuse scattering 
was the chief means of measuring phonon disper- 
sion in solids by fitting lineshapes with expressions 
like Eq. (114). 

3.6 Quasielastic Scattering 

Quasielastic scattering is energy-resolved scatter- 
ing centered on (o = and is the result of neutrons 
interacting with purely dissipative excitations, 
which can be viewed as motions at imaginary 
frequencies. Usually this means scattering by diffus- 
ing nuclei. Typically, quasielastic and inelastic 
scattering are well-separated from each other, since 
vibrational angular frequencies tend to be much 
larger than diffusive jump rates. Sometimes, how- 
ever, analysis of the scattering near w = into its 
elastic, quasielastic and inelastic components is not 
easy. 

A prototype for quasielastic neutron scattering 
is incoherent quasielastic scattering from hydrogen 
in metals. At high temperatures the correlation 



function, G\m(r,t) in Eq. (30), is well-approxi- 
mated by the classical limit, 

GSlr(r,0=Pseif(r,0 

= {8{r+Ro(0)-Ro{t))) , (117) 

which is the finite-^ generalization of Eq. (58). 
Application of jump diffusion theory to the calcula- 
tion of Px\f{r,t) then leads to the Lorentzian 
quasielastic lineshape. 



Sinc(Q,(0) = 






(118) 



where r{Q) depends on the diffusion constant and 
on the structure visited by the diffusing particle. In- 
coherent quasielastic scattering coexists with inco- 
herent elastic scattering given by Eq. (57), which is 
called the elastic incoherent structure factor in the 
concerned literature. This "EISF" is proportional 
to the reciprocal of OxH, the volume accessible to 
the diffusing particle. Thus, the less localized the 
self-diffusion, the smaller the weight of incoherent 
elastic scattering. 

Coherent quasielastic scattering, which was 
alluded to in Sec. 3.3, is important in fluid systems, 
indeed wherever diffusive transport of coherently 
scattering nuclei occurs. Quasielastic scattering also 
is produced by diffusive rotational modes in plastic 
crystals and caged systems. See [4] for an extensive 
treatment of quasielastic scattering. 

4. Neutron Refraction and Reflection 

The Born approximation neglects interference 
between the incident and scattered waves. This usu- 
ally is justified in standard beam-target geometries 
for scattering at wide angles from small samples. 
These interactions become important, however, as 
scattering becomes concentrated into the forward 
direction. Then the sharp distinction between inci- 
dent and scattered beams within the scattering 
medium is lost, and it is necessary to solve the 
Schrodinger equation ("dynamical" scattering) 
without resort to the Born approximation ("kine- 
matical" scattering). 

An important problem where dynamical scatter- 
ing theory is required is a neutron beam incident on 
a semi-infinite scattering medium having a smooth, 
planar boundary, say at z=0. The Schrodinger 
equation for this case is 

^VV(r) + F(r)^(r)=£o^(r), (119) 
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where Eo = A^ko/2m , and the potential is the step 
function, 



V(r) = 



0, 



<I^>cell , 



z > 
z € 0- 



(120) 



The amplitude of the reflected wave works out to 



Mr = 



sin^— '\/sin^(f)—sin^(f)c 
sin^+ \/sin^ (f> - sin^ (f>c 



(126) 



The effective scattering potential in the medium is 
the average of V(r) over the crystalline unit cell. 
This is a good approximation for neutrons beyond 
the Bragg cutoff. In detail. 



(K)cell = P (icoh)cell , 



(121) 



where p =N/V, and 6coh is arithmetically aver- 
aged over the cell. The trial solution is 



'/'('•) = 



r e'*»''+M,e'*'''' z>0 



where ^c is the critical angle, 



(f>c = arcsin \/l— n' 



(127) 



For (f>^(f>c, Eq. (126) shows that |Mr| = l, which 
indicates that the incident wave is totally reflected. 
In this case the normal component of km is imagi- 
nary, and the wave in the medium is evanescent. 
When (f) > (f>c, km is real, and its direction is given 
by Snell's law. The index of refraction in Eq. (125) 
can be written as 



1 Ao 
"-^-2^-n ' 



(128) 



where k[=ko. The subscript "r" denotes the 
reflected beam, while "m" refers to the wave within 
the scattering medium. The angle of incidence ^ 
and the angles of the reflected and refracted waves 
({)[ and ^m, respectively, are conventionally mea- 
sured from the flat in neutron optics. Since the 
model potential is everywhere independent of j: and 
y, the wavevectors ko, k; and km have equal compo- 
nents parallel to the surface. Thus the reflection is 
specular, 0r = . The boundary conditions at z = 
are 



l + Mr = "m , and 
(1 —Ut)ko sin ({) = Mra^ra sin (f)„ 



(123) 



It follows easily from Eq. (122) in Eq. (119) that 



where rj =p(Z>coh)ceii is the average scattering length 
density of the medium and generally is positive. 
Typically n =1—0(10"') for thermal neutrons, 
which gives critical angles 



^c 



^ IT 



(129) 



of only fractions of a degree. The longer wave- 
lengths of cold neutrons give substantially larger 
critical angles, which are required for practical 
applications of neutron guides. For media in which 
t; <0, the index of refraction is greater than unity, 
(f)c is imaginary, and the incident wave is never 
totally reflected (except at ^ = 0). 

The scattering wavevector for neutron reflection 
is defined as 



Q =ko-kc 



(130) 



Pki fi^kl 



2m 



2m 



+ <»^>c 



(124) 



so that Q=2kQsin ({>. The neutron reflectivity, 
R =\u[\^, then is obtained from Eq. (126) as 



This gives the index of refraction of the scattering 
medium, 



R(Q) = 



q^VqEqI 



(131) 



"-Y.-V^- 



{n 



cell 



(125) 



where Qc = 2kosin (f)c. This is equivalent to the 
Fresnel formula for reflectivity in optics. Thus, 
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R{Q) = l for Q^Qc, while for Q>Qc, the reflec- 
tivity falls as 



R{Q) 






(132) 



The asymptotic behavior in Eq. (132) is identical to 
that derived from the Born approximation for the 
potential in Eq. (120) when da/dH, Eq. (60), is 
converted to the reflectivity, an observation made 
only recently by Sinha, et al. [16]. The reflection 
laws from heterogeneous and nonuniform surfaces 
and films depart significantly from Eq. (131) and 
are sensitive to model details. For expositions of 
neutron reflectometry and optics, see [9], [12], 
and [15]. 



[15] V. F. Sears, Neutron Optics, Oxford University Press, New 

York (1989). 
[16] S. K. Sinha, E. B. Sirota, S. Garoff, and H. B. Stanley, 

X-ray and Neutron Scattering from Rough Surfaces, Phys. 

Rev. B38, 2297 (1988). 
[17] K. Skold and D. L. Price, eds., Methods of Experimental 

Physics, Volume 23, Neutron Scattering Academic Press, 

Orlando (1986). 
[18] G. L. Squires, Introduction to the Theory of Thermal 

Neutron Scattering, Cambridge University Press, 

Cambridge (1978). 

About the author: Norman F. Berk is a theoretical 
physicist in the Reactor Radiation Division of the 
NIST Materials Science and Engineering Laboratory. 
The National Institute of Standards and Technology 
is an agency of the Technology Administration, U.S. 
Department of Commerce. 



5. References 



[1] G. E. Bacon, Neutron Diffraction, 3rd ed., Oxford Univer- 
sity Press, Oxford (1975). 

[2] G. E. Bacon, Neutron Scattering in Chemistry, Butterworths, 
London and Boston (1977). 

[3] E. Balcar, and S W. Lovesey, Theory of Magnetic Neutron 
and Photon Scattering, Oxford University Press, New York 
(1989). 

[4] M. B6e, Quasielastic Neutron Scattering, Adam Hilger, 
Bristol and Philadelphia (1988). 

[5] P. A. Egelstaff, ed.. Thermal Neutron Scattering, Academic 
Press, New York (1965). 

[6] G. Kostorz, ed.. Neutron Scattering, Treatise on Materials 
Science and Technology, Volume 15, Academic Press, New 
York (1979). 

[7] M. A. Krivoglaz, Theory of X-Ray and Thermal-Neutron 
Scattering by Real Crystals, Plenum Press, New York 
(1969). 

[8] W. Marshall and S. W. Lovesey, Theory of Thermal 
Neutron Scattering, Clarendon Press, Oxford (1971). 

[9] J. Penfoid and R. K. Thomas, The Application of the 
Specular Reflection of Neutrons to the Study of Surfaces 
and Interfaces, J. Phys. Cond. Matt. 2, 1369 (1990). 
[10] D. L. Price and K. Skold, Introduction to Neutron Scatter- 
ing, in Methods of Experimental Physics, Volume 23, 
Neutron Scattering, K. Skold and D. L. Price, eds.. 
Academic Press, Orlando (1986) p. 1. 
[11] E. Prince, Mathematical Techniques in Crystallography 
and Materials Science, Springer- Verlag, New York (1982). 
[12] T. P. Russell, X-ray and Neutron Reflectivity for the Inves- 
tigation of Polymers, Materials Science Reports 5 (1990) 
p. 171. 
[13] R. Scherm, Fundamentals of Neutron Scattering by 

Condensed Matter, Ann. Phys. 7, 349 (1972). 
[14] V. F. Sears, Neutron Scattering Lengths and Cross 
Sections, in Methods of Experimental Physics, Volume 23, 
Neutron Scattering, K. Skold and D. L. Price, eds. 
Academic Press, Orlando (1986) p. 521. 



30 



